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New Hamiltonian constraint operator for loop quantum gravity 
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A new symmetric Hamiltonian constraint operator is proposed for loop quantum gravity, which is well defined 
in the Hilbert space of diffeomorphism invariant states up to non-planar vertices with valence higher than three. 
It inherits the advantage of the original regularization method, so that its regulated version in the kinematical 
Hilbert space is diffeomorphism covariant and creates new vertices to the spin networks. The quantum algebra 
of this Hamiltonian is anomaly-free on shell, and there is less ambiguity in its construction in comparison with 
the original method. The regularization procedure for this Hamiltonian constraint operator can also be applied to 
the symmetric model of loop quantum cosmology, which leads to a new quantum dynamics of the cosmological 
model. 

PACS numbers: 04.60.Pp, 04.60.Ds 


The singularity theorem of general relativity (GR) is a 
strong signal that the classical Einstein’s equations can not 
be trusted when the spacetime curvature grows unboundedly. 
It is widely expected that a quantum theory of gravity would 
overcome the singularity problem of classical GR. A very les¬ 
son that one can learn from GR is that the spacetime geometry 
itself becomes dynamical. To carry out this crucial idea raised 
by Einstein 100 years ago, loop quantum gravity (LQG) is 
notable for its nonpertuibative and background-independent 
construction [1^]. The kinematical Hilbert space of LQG 
consists of cylindrical functions over finite graphs embedded 
in the spatial manifold. The quantum geometric operators 
cotTesponding to area [5, 6], volume [5, 7, 8], length[9-l 1], 
ADM energy [12] and quasi-local energy [13], etc. have dis¬ 
crete spectrums. The LQG quantization framework can also 
be generalized to high-dimensional GR [14] and scalar-tensor 
theories of gravity [15, 16]. A crucial topic now in LQG is 
its quantum dynamics, which is being attacked from both the 
canonical LQG and the path integral approach of spin foam 
models. In the canonical approach a suitable regularization 
procedure was first proposed by Thiemann to obtain well- 
defined Hamiltonian constraint operators [17]. The Hamil¬ 
tonian constraint operators obtained in this way will attach 
new arcs (edges) and hence create new trivalent co-planar 
vertices to the graph of the cylindrical function which it acts 
[17, 18]. The quantum dynamics determinatted by the Hamil¬ 
tonian constraint operator is well tested in the symmetric mod¬ 
els of loop quantum cosmology (LQC) [19]. The classical 
big bang singularities are resolved by quantum bounces in the 
models [20-22]. However, there are ambiguities in the graph- 
dependent triangulation construction of this operator. In or¬ 
der to obtain the on shell anomaly-free quantum algebra of 
the Hamiltonian constraint operator [23], one has to employ 
degenerate triangulation at the co-planar vertices of spin net¬ 
works in the regularization procedure of the Hamiltonian.^ 
This treatment implies that the regularization procedure has 
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essentially neglected the Hamiltonian at the co-planar vertices 
before acting the regulated operator on them. Otherwise, this 
kind of Hamiltonian constraint operator would generate an 
anomalous algebra in the full theory, unless one inputs certain 
unnatural requirement to the interaction manner of the edges 
of the graph and the arcs added by the Hamiltonian opera¬ 
tor [24]. The Hamiltonian constraint operators proposed re¬ 
cently in [25-27] do not generate new vertices on the graph of 
the cylindrical function and hence are anomaly-free on shell. 
However this kind of action can not match the quantum dy¬ 
namics of spin foam models where new vertices are unavoid¬ 
able in their construction [28]. A regularization of the Hamil¬ 
tonian constraint compatible with the spinfoam dynamics was 
considered in [29]. However, the resulted Hamiltonian oper¬ 
ator acts non-trivially on the vertices that it created and thus 
has still an anomalous quantum algebra. It is therefore natural 
to ask the question whether one can construct some Hamilto¬ 
nian constraint operator with the following properties: (i) it is 
diffeomorphism covariant, symmetric and anomaly-free; (ii) 
it generates new vertices; (iii) its action on co-planar vertices 
is not neglected by some special regularization procedure, and 
there is no special restriction on the interaction manner of the 
edges of the graph and the arcs added by its action. We will 
show that the answer is affirmative. An alternative quanti¬ 
zation of the Hamiltonian constraint in LQG possessing the 
above three properties will be proposed. The regularization 
procedure of the Hamiltonian operator can also be applied to 
LQC models. 

The Hamiltonian formalism of GR is formulated on a 
4-dimensional manifold M - R x L, with L being a 3- 
dimensional spatial manifold. In connection dynamics, the 
canonical variables on E are the SU{2)- connection A'^ and 
the densitized triad E'j, with the only nontrivial Poisson 
bracket (A]j(x), £*(y)] = Kf]6^{x,y), where k = %nG and yS is 
the Barbero-Immirzi parameter. The Hamiltonian constraint 
reads 
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where = 25[aA]^j + is the curvature of A‘^, K'^ 
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is the extrinsic curvature of 2, and det (q) is the determinate 
of 3-metric q^b = with e‘^ being the co-triad. H^(N) 

and T{N) are called the Euclidean and Lorentzian terms of 
the Hamiltonian constraints respectively. Both of H^(N) and 
T{N) depend on the canonical variables in non-polynomial 
ways. Besides the indication of spin foam models, it is ar¬ 
gued in [30] that the momentum variables in H^(N) also im¬ 
ply the creation of new vertices by its action. Thus we adopt 
the so-called semi-quantized regularization approach devel¬ 
oped in [31] to derive a new Hamiltonian constraint operator, 
which creates new vertices as well. The Hamiltonian is not ne¬ 
glected at the co-planar vertices of spin networks by the regu¬ 
larization. But the result of its action on the co-planar vertices 
is zero. Hence it has an anomaly-free algebra on shell. 

Let us first consider H^{N). By introducing a characteristic 
functionsuch that \\mxe{x,y)l- 6^{x,y) and using 

the point-splitting scheme, it can be regularized as 

//^(iV)=llim r d^xN{x)V7^'hquFUix)E‘‘{x) 

ZK e-^0 Jj; ^ ' 


x£d^yxAx,y)E^b(y)V^f’ ( 2 ) 


where V(x,^) := y/det(q)(x). Since the volume operator has 
a large kernel, the naive inverse volume operator is not well 
defined. However, one can use the idea in [32] to circumvent 
this problem by defining a permissible inverse square root of 
volume operator as 
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(3) 


where is the standard volume operator in LQG (see [7]) 
corresponding the volume of the cube with center y and radial 
e. It is easy to see that qualitatively j-jjg same proper¬ 

ties with y. Thus we can promote the classical volume in (2) 
into its quantum version (3) and replace both of densitized tri¬ 
ads in (2) by corresponding operators E'^iy) - ISA^^^iy) 

where - tiK. Acting on a cylindrical function fy, the result 
formally reads 
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where 0(t,t') - 1 for f' > t and zero otherwise, X'^(t) ti[(he(o,t)Tkhe^,j)Yd/dhe(ox], \- 

^{{he(Q,t')Tjhe{i\i)Tkhe{t,\)y here Tk -^cTk with cTk being the Pauli matrices, and T denotes transpose. Parti¬ 
tioning of the domain [0,1] as A segments by inputing A - 1 points, 0 = fo, fi, ■ ■ ■ , fv = 1, and setting Af„ s - f„_i = S, 
the integral in (4) can be replaced by the Riemann’s sum in a limit. Then (4) reduces to 
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Since the volume operator and hence vanish at divalent 
vertices, for sufficiently small e, the only non-vanishing terms 
of the summation in (5) correspond to those of m - n - 
Moreover, the terms corresponding Xom - n - 1 in the second 
summation of (5), which involves only summation over e{- 
e'), vanish due to e,yi:A^^(e(fQ))e“(fg)e*(fo) = 0. For e + e', we 
have 

~ ^sgn(e',e)trf(/!„^,^T,), (6) 

where Nc - ^ half-integer represent¬ 


ing a spin representation of SU{2), Oe'e is a loop formed 
by adding an arc between e'{6) and e{S), and sgn(e',e) := 
sgn[eafee'‘’(0)e*(0)] is the orientation factor which can be pro¬ 
moted into its quantum operator. From the property of V, we 

know that is independent of e. Thus we can take 

the trivial limit e ^ 0 and obtain 
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where Ny = N{v) and J‘^ = -iX^, is the self-adjoint right- 
invariant operator. The assignment of a^’e is diffeormorphism 
covariant in the sense that, for any cp 6 Diff(2) there exists 
ip' e Diff(E) such that ip'{ip{y)) = ip{y) and ip' = 

ipiae'e) [17]. Applying on the n-valent non-planar ver¬ 
tex V of T'' the intertwiner iy associated to v will be changed 
y,j,i 

to /[, while, as a multiplication operator, ha^,^ will change the 
spins associated to the segments of e' and e. Hence the action 

of ^le'e is by 



where_/ e_ {|/ - ,f + {}, j e {\j + t], and 

//(/, /, j, j, - ■ ■)aie coefficients with “ ■ ■ ■ ” denoting spins as¬ 
sociated to other edges incident at v. 

The above regularization approach can be similarly used 
to quantize the Lorentzian term in the Hamiltonian (1). The 
Lorentzian term can be regularized as 

^(A^) = lim £ A^(^) 

x£d^yxAx,y)E^(y)V^lf. (9) 

Following the regularization procedure of the Euclidean term, 
Eq. (6) is now replaced by 

d^Kt-^(e'(0))K^\e'(0))e'‘‘(0)e\0) ^ sgn(e', e)K^^,K^, ( 10 ) 


where K^, with as the start¬ 

ing segment of e' with parameter length 6, and we have 
used the identity Ki = ^[ALK] = -^txeiTj{Aa, K]), with 
K ^^d?xK'^E" - ^[El^(l), V) [17]. Replacing the Poisson 
brackets by commutators times ll(iH), we obtain the quan¬ 
tized Lorentzian term 


TsiN) ■ fy = 


2(1 +yS2) 


2k 


’-L ^ N.V. 

veV(y) 
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J] sgn(e',e)Kj,KUi,4 


V;'^^-fy, (11) 


Notice that for a given y, there are indeed finite terms con¬ 
tribution to the summation no matter how fine the partition 
is. Hence the operator (13) is well defined in Since the 

volume operator and hence vanish at internal gauge in¬ 
variant trivalent vertices as well, Hs{N) acts nontrivially only 
on non-planar vertices with valence higher than three. Note 
that we can also employ a fourfold point-splitting scheme as 
in [31] to regulate the Hamiltonian constraint. In this case 
1/ -^det(^) in (1) can be absorbed into Poisson brackets, and 
the final operator takes the same form as (7) and (11), except 

that Vy is replaced by 


y-l/2 


4x8 

3!£(v) 



J;ns/nj*;=V 


(14) 


where It can be shown that quali- 

tatively has the same key properties with Vy^^^ [33]. 

Let <t) Cyl“(/7[/^) be the set of smooth cylindrical func¬ 

tion and ‘Fpjg: be the space of diffeomorphism-invariant distri¬ 
butions on <I). Then the number ip{Hs{N) ■ /) for any / 6 (D 
and ^ 6 fhQjjf depends only on the diffeomorphism class of 
the loop assignments. Hence the limit d —> 0 can be taken in 
the Rovelli-Smolin topology as [24, 34] 

[{H{N))' ■ <p)(fy) := lim((H,(/V))' ■ (A)(/y) 

:^\imil/(Hs(N)-fy). (15) 

However, {H{N)y is not well defined in since it does 

not keep Opjjf invariant. Taking account of the fact that the as¬ 
signment of loops by H^iN) is diffeomorphism covariant and 
it acts nontrivially only on non-planar vertices with valence 
higher than three, we consider almost diffeomorphism invari¬ 
ant states which are obtained from the spin network states by 
averaging over their images under diffeomorphisms but leav¬ 
ing fixed sets of non-planar vertices with valence higher than 
3 in the spatial manifold invariant, parallel to the proposal in 
[25]. Given a graph y, we denote its non-planar vertices with 
valence higher than 3 by ynp 4 (y), the group of all diffeomor¬ 
phisms preserving ynp 4 (y) by Diff(E)v_p^(y), and the diffeomor¬ 
phism acting trivially on y by TDiff(S)^. Eor any fy 6 O, we 
define a map 77 ; <1) ^ O' by 

Yj V^-fy, (16) 


where 

/7;/,l]),withl:=^[H,"(i),y]. 

( 12 ) 

Hence the total regulated Hamiltonian constraint operator is 
given by 

mm ■ /y ;= ^ m [H^ - t) • /^ =; J] Nyfly ■ fy. 

veV(y) veV(y) 

( 13 ) 


where is the unitary representation of <p, and Ny is a nor¬ 
malization factor. We can equip the space Onp 4 7/(0) with a 
natural inner product as (77(/)|77(g)) ;= T]{f)(g),'^f,g e O. The 
new Hilbert space 'Knp 4 is defined as the completion of Onp 4 . 
Note that Onp 4 is also in the dual space of the space of diffeo¬ 
morphism invariant states up to all vertices of spin networks 
which was fully discussed in [25]. Whether (HiN))' can be 
well defined in '77np4 is a delicate issue. In particular, if the 
action of H^iN) on a non-planar vertex with higher valence 
could deduce the valence so that it becomes less than 4 or the 
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vertex becomes co-planar, (H(N)y would be ill-defined on the 
resulted almost diffeomorphism invariant states.^ Fortunately, 
we can use the freedom of choosing the spin representations 
{ attached to each new added loop in ( 6 ) to ensure that the 
valence of any vertex would not be changed by the action of 
HsiN). Then it is straightforward to see that (H{N)y is well 
defined in 'Knp 4 . By the dual action, (H(N)y will annihilate 
the arcs like the one connecting V 2 and V 3 in ( 8 ). Since Hg(N) 
vanishes at co-planar vertices, it is obvious that 

{[(H(M)y,(H(N)y]-4>)(fy) 

^cj, {{hANWM) - HAM)Hs(N)) ■ fy) 

=0, V/^e(l),0e<h„p4. (17) 

Hence for any ^ e which is also a distribution on Onp 4 , 
we have 

(([(H(M))', {HiN))']) ■ -a) (<A) = <A([(77(M))', (H(M))')] ■ <a) 

= 0, V,A e <l>„p4. (18) 


Therefore the quantum algebra of the Hamiltonian constraint 
operators {H{N)y is anomaly-free on shell. A symmetric 
Hamiltonian constraint operator corresponding to {H{N)y can 
be defined in ' 7 /np 4 as 

(^sym(^)), _ 1 ((H(A))' + ((H(N)yy ) , (19) 

where ((//(A))')^ denotes the adjoint of (HiN))'. The action 
of ((H(A))') on an almost diffeomorphism invariant state will 
create co-planar vertices and arcs like in ( 8 ). Thus it is easy to 
see that the quantum algebra of is also anomaly- 

free. 

Now we test the above regularization technique of 
the anomaly-free Hamiltonian constraint operators in the 
symmetry-reduced model of LQC. For simplification, we con¬ 
sider only the spatially fiat and isotropic model. The con¬ 
struction can be generalized to the other cosmological mod¬ 
els directly. Introducing an elementary cell 'V adapted to the 
fiducial triad, the connections and the densitized triads Ef 
can be expressed as A'^ - °e‘^ and Ef - 

where {°ef, ’’ef) are a set of orthonormal co-triads and triads 
compatible with the fiducial metric [19]. The basic (non¬ 
vanishing) Poisson bracket is given by jc, p) = y. To pass 
to the quantum theory, one constructs a kinematical Hilbert 
space TY®™ = L^(RBohr, dpBohr), where RboHi is the Bohr com- 
pactification of the real line and dpBohr is the Haar measure on 
it [19]. The holonomy of Af along an edge, parallel to the 
triad vector °e?, of length with respect to °q^i, is given 

by hf - cos + 2 sin where I is the identity 2x2 ma¬ 
trix. Because of homogeneity, the Hamiltonian constraint can 
be written as [19, 20 ] 


1 F 1 

H = —=-T 

/32 2 /c/?2 


X‘ 


f P' P" pb 
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d^X 


Vdet(q') 
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Since -\Jdet{q) - v(Vg^ -sjdsti^q)^ where V is the 

physical volume of 'V, we can write H as 

Note that Ef can be directly quantized as 

= y;2/3 ^detC’q) ’’efp ;= yJdet(oq) ‘>ef x (-inf) ^ ■ 

( 22 ) 


After the action of Ef and Ef, we can integrate over the 
plane spanned by °ef and °ef to yield a holonomy along a loop. 
Promoting functions into corresponding operators, the regu¬ 
lated operator corresponding to ( 20 ) reads 




1 

2k0^N{ 


_ eijkykt\hg/‘Ti\ 

y-1/2- \^:Elppv-ii 2 ^ 


(23) 


where oAj, is a loop along a square in the j-k plane spanned by 
a face of 'V, each of whose sides has length with respect 
to "qab, and = hfhf . Since the area operator 

in LQG has a minimum nonzero eigenvalue 2 which 

introduces a nature cut off to the size of the loop we obtain 

pP\p\^2^nl3£l. (24) 

The symmetric Hamiltonian constraint operator can be de¬ 
fined by 


^.sym _ l(^^^t) (25) 

This Hamiltonian constraint operator is obviously simpler 
than those appeared in LQC models (e.g. the APS Hamil¬ 
tonian operator [20]). It is straightforward to check that the 
Hamiltonian operator (25) has correct classical limit and the 
classical big bang singularity can also be avoided by a quan¬ 
tum bounce of this dynamics. 

To summarize, the Hamiltonian constraint of GR is suc¬ 
cessfully quantized in LQG by adopting the semi-quantized 
regularization approach. The resulted operator is symmetric 
and well defined in the Hilbert space ' 7 fnp 4 of diffeomorphism 
invariant states up to non-planar vertices with valences higher 
than 3. The action of this Hamiltonian constraint operator cre¬ 
ates new trivalent co-planar vertices to the spin networks but 
does not change the valence of the acted vertices, and hence 
it can be symmetric and match the quantum dynamics of spin 
foam models. Meanwhile, our regularization procedure en¬ 
sures that there is less ambiguity in the construction of the 
operator. The Hamiltonian is not neglected at the co-planar 
vertices of spin networks by the regularization. But the re¬ 
sult of its action on the co-planar vertices is zero. Hence the 
quantum algebra of the new Hamiltonian constraint operator 
is anomaly-free on shell. It is thus possible to find solutions 
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of the quantum Hamiltonian constraint in 'Knp 4 . The regu¬ 
larization procedure for the Hamiltonian operator in full the¬ 
ory is also applied to the isotropic model of LQC. The re¬ 
sulted Hamiltonian constraint operator is simpler than those 
appeared in LQC models. Meanwhile, it inherits the qualita¬ 
tive properties of the previous Hamiltonian operators in LQC, 
so that it has correct classical limit and the classical big bang 
singularity can also be avoided by a quantum bounce. 
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